Effect of Electric Field on One-Dimensional Insulators: A DMRG study 
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We perform density matrix renormalization group (DMRG) calculations extensively on one dimen- 
sional Mott and Peierls chains with explicit inclusion of the static bias to study the insulator-metal 
transition in those systems. We find that the electric field induces a number of insulator-metal 
transitions for finite size systems and at the thermodynamic limit, the insulating system breaks 
OO ■ down into a completely conducting state at a critical value of bias which depends strongly on the 

' insulating parameters. Our results indicate that the breakdown, in both the Peierls and Mott in- 

' sulators, at thermodynamic limit, does not follow the Landau-Zener mechanism. Calculations on 

various size systems indicate that an increase in the system size decreases the threshold bias as well 
f— I ' as the charge gap at that bias, making the insulator-metal transition sharper in both cases. 
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Strongly correlated low-dimensional electronic systems have attracted much interest because of their unique low- 
^ , energy characteristics[l|,[2|, [^[3, Q. These systems are almost alvifays insulators due to various interactions in reduced 
I ■ dimensions and are commonly described by Hubbard, Peierls or related Hamiltonians[l, While the ground state 
|H ' of a half-filled Hubbard system is a Mott insulator, the electron-lattice interactions lead to Peierls instability in 
^ ■ low-dimcnsional systems. 

The effect of electric field on such systems has generated much interest in recent years due to the practical interest 
in tuning their dielectric and piezoelectric properties, more so for the Mott insulators. Many electronic conduc- 
tion processes seem to suggest electric field induced phenomena, such as negative differential resistance in molec- 
ular electronics B [^. Recent experiments on low-dimensional Mott insulators, namely, Sr2Cu03, SrCu09 [lOt and 
• La2-2:Sr2:Ni04[ll| suggest a dielectric breakdown in presence of an external electric field. While the Peierls case is 
Q ' relatively easier, a tractable computational method which takes into account the static electric field and its response 
on the excited correlated electronic states is still lacking. One of the earliest theoretical approaches considering electric 
field was the Bcthe-awsate methodp^.[l3|. where, asymmetric tunneling terms towards the left and right were consid- 
ered inordcr to describe the non-equilibrium situation and hence capture the flavor of an electric field. However, the 
^ ' imaginary gauge term giving rise to a non-hermitian Hamiltonian makes it very difficult to relate this to a system in 
presence of a real electric field. In order to account for this, another approach to this problem has been to employ the 
time-dependent Schrodinger equation to study the time evolution of the many-body wave functions and levels, where 
2. ] the electric field is applied via a time-dependent Aharonov- Bohm fluxQ." However, in this case the calculations 
are restricted to only ring structures (to avoid ambiguities arising from the electrodes) and for system sizes < 10 
sites. Recently, however, there have been a number of studies using the real time evolution of the ground state in the 
presence of a source drain potential [l5|. However, a microscopic understanding of the insulator-metal transition and 
the quantification of the critical static electric field, required to induce it in real extended systems is still lacking. 

For a clear understanding of the effect of electric field on finite size systems as well as in the thermodynamic limit, 
in this letter, we use the Density Matrix Renormalization Group (DMROflGl. [iTj method which is known to be 
highly accurate for low-dimensional interacting systems. For the first time, we have included static electric field in the 
DMRG algorithm and have obtained ground and low-energy eigen states of a one-dimensional system with various 
d ] sizes. The static electric field is included as a ramp potential and we find that the electric field can induce a number 
of insulator-metal transitions for finite size systems, with the transition depending on the Hamiltonian parameters. 
We analyze our results at thermodynamic limit in the light of Landau-Zener mechanism [isj . for both the Peierls and 
Mott insulating systems. To get the thermodynamic behavior we consider the total bias instead of field. 

We consider a one dimensional chain representing a conjugated polymer or any other one-dimensional system, 
described by the Peierls-Hubbard Hamiltonian, 

i 

+ ^ X!"*T'^ii (1) 



where t is the hopping term, U is the Hubbard term and i5 is the bond alternation parameter. We set t = 1 and for 
the Mott-insulator, we consider 5 = and nonzero U , while for the Peierls-insulator. U is set to zero with nonzero 
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FIG. 1: {a.)N — > oo values of /i+(square) and /i_ (circle) at zero bias as a function of U. The dotted lines correspond to exact 
Bethe-ansatz results. (b)Converged ground state energy per site for half-filled system (N electrons) vs applied bias (V) for 
U = (X), 1 (circle), 2 (square), 3 (diamond), 4 (triangle) and 5 (star), (c) shows the energy of the systems with N-fl (up 
triangle (square)), and N-1 (down triangle(circle)) electrons for U = 3(5) for a chain of 40 sites as a function of bias. The 
dashed lines represent the corresponding ground state energy (E(N)). 

6. The external electric field applied on the system has the form of a ramp potential, distributed over all the sites 
in such a way that the potential Vi at site i becomes — -j + *-/vTT' ^ticre V is the applied voltage and N is the total 
number of sites in the ID chain. This form of the potential ensures that the bias varies between —V/2 to V/2 across 
the molecule. The potential adds an extra term ^ Via^tai to the above Hamiltonian. 

i 

Our DMRG results compare fairly well with the V ~ Bethe- a^sate ground state energy and charge gap, for all 
values of U , with a density-matrix cut-off m = 140. For nonzero bias, the DMRG results compare up to numerical 
accuracies with the exact diagonalization results of finite sizes up to TV = 16. For U = 0, the problem can be exactly 
solved and we obtain the ground state and excitation spectrum in presence of bias using tight-binding one-electron 
formalism. In the Mott (Peierls) case, U{5) has been varied from to 5 (0 to 1) with bias from to 6 volts in steps 
of 0.1 volt. 

For many-body models, the charge excitation gap is defined as the difference between the energy required to add 
{fj.+ ) and remove (/i-) electrons from the ground state [l^. 

'^charge = M+ ^ (2) 

where ^+ = E{N + 1) - E{N) and /i_ = E{N) - E{N - 1). E{N), E{N + 1) and E{N - 1) are the energies of the 
half-filled system and the systems with one extra and one less electron respectively. 

To understand the energy cost due to the addition or removal of electrons in the presence of bias, we have calculated 
/i+ and /i_ numerically using DMRG for a range of bias. In Fig. 1(a), we have plotted /i+ and ^_ as a function of U 
at zero bias. The finite size DMRG results {fi+ and //_) are extrapolated to the thermodynamic limit {N — s- oo) for 
every U value. For a clear demonstration of the accuracy of our results, we have also shown the Bethe- ansate /Lt+ and 
/i_ values derived for every U in the same plot. As can be seen, our numerical extrapolated values compare fairly 
well with the exact thermodynamic Bethe-a77.sate results. Both ^_|_ and fi- increase with increase of U. However, the 
former increases with larger slope than the latter, resulting in increase of the charge gap with increasing Hubbard 
repulsion. This is because, a higher value of Hubbard repulsion localizes the electrons more and leads to higher charge 
gap. 

The ground state of the Hamiltonian with nonzero ?7 is a spin-density wave insulator with one electron at every 
sitejiOl- In Fig. 1(b) we have plotted the converged (extrapolated value at iV — > oo limit) ground state energy of the 
half-filled system as a function of applied bias for different U values. As can be seen, for a fixed U value, the ground 
state energy remains almost constant upto a certain bias, and beyond it, the system starts stabilizing rapidly with 
applied bias[2l|. It can thus be interpreted from Fig. 1(b) that the applied electric field stabilizes the system after 
a threshold value, required to overcome the effect of the electron repulsion, U. Thereafter, the system gains kinetic 
energy as the electrons start hopping in the direction of bias. For U = 0, the stabilization starts as soon as the bias 
is turned on. 

To understand the response of the excited states of the system with electric field, for a given U, we have computed 
E{N + 1) and E{N — 1). Both the energies for a particular system size {N = 40) are plotted in the Fig. 1(c) as a 
function of bias for two representative U values, together with E{N). As like the ground state, excitation energies 
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FIG. 2; Charge gap vs bias for A'^ = 20(solid line), 26(dashed line) and 40(dotted line) for [7 = 4. Inset shows the charge 
densities on first (solid lines) and last site (dashed lines) for the N-l(a), N(b) and N+l(c) electron systems as a function of 
bias for U — 4(circle) and 5(star) w 
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FIG. 3: Bias at first (circle), second (square) and third (triangle) minimum as a function of for U — 4. Dashed lines show 
the extrapolation. Inset shows the /i+(star) and //-(diamond) at first minimum as a function of 1/N for U = 3(solid lines) and 
4(dashed lines). Dotted lines are extrapolation. Error bars are of the same size as the symbols. 



(both E{N + 1) and E{N — 1)) decrease with increase in bias, however, their difference always remains equal to U 
for aU values of bias, as expected. For a given value of U, although the slope of E{N + 1) and E{N — 1) arc same, 
the slope of E{N) is different leading to the interesting phenomenon of insulator-metal transition, as described later. 
Fig. 1(c) also shows that the stabilization of the system with one extra/less electron needs higher bias for higher value 
of [/. 

To characterise the bias at which this insulator-metal transition occurs accurately, we plot in Fig. 2 the charge gap as 
a function of bias for different finite size systems for a representative value of U. It can be clearly seen that the charge 
gap shows an oscillation with bias, going through a number of minima and maxima. To understand the underlying 
physics, we have calculated the onsite charge densities as a function of bias for various system sizes with several U 
values. We present this in the inset of Fig. 2 for the half-filled state and the states with one extra and one less electron 
than half-filling for two representative values of U for a finite chain with N = 30. For clarity, wc plot charge densities 
at only the first and the last sites. At zero bias, the ground state charge density at every site of the system is the same 
and it remains almost the same with increase in bias upto the bias corresponding to the first minimum of the charge 
gap (Ai?). However, after that, they show a large shift in the direction of bias, giving rise to charge inhomogeneties. 
A higher value of U requires a higher bias to shift the charge density as can be seen from the inset of Fig. 2 and thus 
the bias corresponding to the first AE minimum increases with increasing U. Interestingly, the repetetive period of 
charge gap going through minima with increase in bias is due to the role of charge stiffness. The external bias tends to 
shift the charge densities towards one electrode with the nullification of U at the first AE minimum. However, beyond 
this, an increase in bias results in further hopping of charges leading to double occupancy of more sites, with electron 
repulsion overwhelming the kinetic stablization, thereby increasing the energy gap. Further increase of bias nullifies 
this effective repulsion, resulting in the next charge gap minimum. Hence, such a variation in charge gap resulting in 
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FIG. 4: Vc as a function of U from our results (circle), and Bethe-ansatz prediction of Landau-Zener tunneling (dashed line). 
Inset shows Vc as a function of S. 

near-metallic (charge gap is not zero) behavior in various bias regions is due to the interplay of the Hubbard repulsion, 
finite system size and the spatial gradient of the external bias. Interestingly, as can be seen from Fig. 2, with increase 
in the size, the magnitude of bias corresponding to the first AE minimum reduces and the periodicity of occurence of 
succesive minima thereafter also narrows down. The field (V/N) required to reach the first AE minimum decreases 
more sharply with increase in system size. 

To explore the size dependence of this insulator-metal transition, in Fig. 3 we plot the bias values corresponding 
to the first, second and third AE minimum separately as a function of 1/N and extrapolate the plots to iV ^ cxo 
limit. All the minima converge to the same point at the thermodynamic limit, clearly indicating that at A'^ — > oo 
limit, the system goes to conducting state at some bias which we term as critical bias, Vc- At this Vc, the ground 
state of the system changes its slope and starts stabilizing rapidly as can be seen from Fig. 1(a). In the inset of Fig. 3, 
we have shown and /i_ at first minimum as a function of inverse system size for [/ = 3 and 4. With increase 
in system size, /i+ and /i_ come close to each other and thus the gap at first minimum becomes lower with larger 
deconfinement. From our finite size data analysis, with extrapolation to A^ — > oo limit, we find that with increase in 
Hubbard repulsion term, U, the difference between /i+ and /i_ at first minimum for a particular finite chain decreases, 
as can be seen from the inset. Thus with increase in [/, the length of the chain at which the charge gap vanishes, 
decreases. Hence, as is evident from the inset, for a given value of U, our method enables us to predict the length of 
the chain at which the system will become completely conducting. 

It is clear from the above discussion that, the bias required for the closing of Mott gap depends sensitively on 
the value of U. To understand how Vc depends on U, we plot in Fig. 4, the Vc as a function of U obtained from 
our calculations, along with that calculated using the Landau-Zener formula [2^ on the initial charge gaps obtained 
using the Bethe- ansatz. Fig. 4 shows non- linear dependence of Vc on U. We have also analyzed the insulator- metal 
transition in case of Peierls insulators. In this case, the difference between the highest occupied and lowest unoccupied 
single electron states (opposite parity) is the charge gap. Here too, the critical bias has a sensitive dependence on the 
bond-alternation parameter, 6, as has been described previously by us[2^, H^. Inset of Fig. 4 displays the variation of 
converged Vc with 6. To now compare and contrast the critical bias obtained for the Peierls and Mott insulators in light 
of the Landau-Zener formula, we analyze the initial charge gap and the Vc dependence on the model parameters, S and 
U. For Peierls insulators, the V = charge gap is linear in 6, and the critical bias Vc, calculated at the thermodynamic 
limit (inset in Fig. 4) has a near linear dependence on S. This indicates that the Landau-Zener formula which has 
the form Vc{S) oc [Achargei^)]'^ , is not obeyed although previous studies on finite size systems indicated that tunneling 
occurs between the Landau quasi-degenerate valence and conduction levels [H, [l^, [2^, resulting in an exchange of 
their symmetries. This is because, at the large A^ limit, the gap at the critical bias goes to zero and the system goes to 
a completely conducting state. In case of the Mott insulator, the initial charge gap has a nonlinear dependence on U, 
and as seen in Fig. 4, the critical bias also shows a nonlinear dependence on U. A fitting gives Vc{U) oc [Acharge{U)]" , 
with n ^ 1, which, unlike the results from previous studies |22| clearly indicates that, the tunneling between many 
body levels does not follow Landau-Zener mechanism in thermodynamic limit. Note that, in our case Vc corresponds 
to the total bias applied to the chain with A ^ oo, while the time dependent studies [l^, [2^ calculate the critical 
field (Vc/N) by converging it upto a certain finite chain length. We do not consider the effect of polarization on 
the applied electric field in our calculations. Because the inclusion of polarization does not change the nature of the 
insulator-metal transition as we have observed earlier [23. It can only change the quantitaive estimation of the bias 
at charge gap minima to some extent. Moreover, in presence of Hubbard repulsion term U the ramp nature of the 
electric field is retained even with inclusion of polarization effects in our calculations p6|. 
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Qualitatively, the slope change of the ground state can be visualized even for a 2-sites Mott-insulating system. 
For V = 0, the ground state of the half-filled system comprising of four basis states, has large contibution from the 
singly occupied sites, with second order [t^ /U) contribution from the double occupancy sites. However, the external 
bias nullifies the Hubbard repulsion, leading to a mixing of the double and single occupancy sites. The ground state 
charge density shifts in the direction of the applied electric field, aXV^ 3U. On the other hand, the state with one 
extra electron comprising of two basis states, stabilizes with a slope which docs not depend on U. This is because the 
external bias results in the hopping of electron from one site to another, leading to a preference of one of the basis 
states as compared to the other. But since both the basis states are energetically degenerate, this hopping occurs 
without any Hubbard energy cost, rather the energy reduces due to kinetic stabilization. This results in a slope that is 
independent of the Hubbard repulsion. Additionally, the Hubbard repulsion also does not play any role on the 2-sites 
one electron case. The 2-sites problem thus captures the ground state behavior, and not the transition pheonomena, 
since the interesting features described earlier [AE oscillation) leading to breakdown behavior occur primarily due to 
the response of the excited states to an applied static electric field. 

In conclusion, we have shown that, in thermodynamic limit, the insulator-metal transition in presence of an external 
static electric field does not follow Landau-Zener mechanism in both Peierls and Mott insulators. Our DMRG 
calculation is the first of its kind, which allows us to study various system sizes with better understanding of the 
general insulator-metal transition, with quantitative predictability. 
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CSIR, Govt, of India. 
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